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, We consider chiral symmetry breaking at nonzero chemical potential and discuss the 

f**^ ■ relation with the spectrum of the Dirac operator. We solve the so called Silver Blaze 

"**>^ ' Problem that the chiral condensate at zero temperature does not depend on the chemical 

potential while this is not the case for the Dirac spectrum and the weight of the partition 

function. 
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1. Introduction 



The QCD phase diagram has been explored in great detail during the past decade. 
' New phases have been discovered and a variety of parameters have been considered 

| 1 | ^ I 3 | 4 | ! | > |6JZ| i n particular, a great deal of progress has been made for QCD at large 
nonzero baryon chemical potential (denoted by /i) . Based on model calculations and 
asymptotic expansions, it is now generally accepted that the ground state of QCD 
in this region is in a phase with a diquark condensate in a color-flavor locked stated 
The QCD phase diagram is best known at /i = where it has been explored exten- 
sively by lattice simulation^. At a critical temperature of about 160 MeV we expect 
a cross-over from a phase with a nonzero chiral condensate to phase with restored 
chiral symmetry. At zero temperature, the QCD partition function is independent 
of until it has reached the lightest excitation with nonzero baryon number. This 
means that the chiral condensate remains constant up to /i = /3 — e (with e the 
nuclear binding energy). The transition to the chirally restored phas e is believed to 
be of first order at low temperatures and ends in a critical endpoinlPn^^EEIl A 
schematic QCD phase diagram is given in Fig. 1. In this lecture we will consider 
the QCD partition function at zero temperature for chemical potentials [i <C tojv/3, 
where QCD can be described by a chiral Lagrangian. 
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Fig. 1. Schematic phase diagram of the QCD partition function at nonzero temperature and 
chemical potential. The short curve denotes the nuclear matter liquid gas transition. 



2. The QCD Partition Function 

The QCD partition function at temperature 1/(3 is given by 

Z QCD ( M ,T)=^ e -^-^), (1) 

k 

where the sum is over all quantum states of QCD with energy Ek and quark number 
iVfc. The partition function is a monotonously nondecreasing function of the chemical 
potential. For T = its value remains constant upto/i = TOjv/3— e (see Fig.|2J). This 
implies that thermodynamic observables, such as for example the chiral condensate, 
do not depend on the chemical potential up to this point. However, this simple fact 
becomes a mystery when we consider the Euclidean QCD partition function given 
by (with quark masses denoted by m/) 

N f N f 

ZqcdG", T) = <n det(D + M7o + m f )) = ( jj []( Afe + ^ 

/=! /=! k 

Here and below, the brackets denote the average over gauge field configurations 
weighted by the Euclidean Yang-Mills action. The quark determinant has been 
expressed in terms of the eigenvalues of the Dirac operator given by 

(D + ^7o)0fe = A fc 0fc. (3) 
The chiral condensate is thus given by (we only consider m/ = m from now on) 

1^)1 = lim^^logZQCD^T) = limi / £ ^-L-A (4) 

* \ k k I QCD 

The mystery is that the eigenvalues and the quark determinant strongly depend on 
the chemical p oten tial but the partition function and the chiral condensate do not. 
Thomas CoheifEl coined this problem as the Silver Blaze problem after the title of 
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Fig. 2. The QCD partition function as a function of the chemical potential for zero and nonzero 
temperature in the absence of nuclear binding energy. The temperature of the critical endpoint in 
Fig. 1 is denoted by T en( j. 



a novel of Conrad Doyle where a racing horse with the name of the title disappeared 
without the dogs barking. The aim of the lecture is to elucidate this problem. 

3. Solution of the Silver Blaze Problem 

The spectral "density" of the Dirac operator is given by 

p(z,z*,m,p) = (J2s 2 (z-Xk)det(D + ^ + m) N f). (5) 

k 

The quenched spectral density is obtained for Nf = 0. We also use the resolvent 
defined by (here and below we take the volume V finite but large) 

G(z, z*,m, p) = (1 V — L- det(D + p 7o + m) N f), (6) 

k 

so that (notice that p{—z, —z*, m, p) — p(z, z* , to, p)) 

V 7T 

Since the chiral condensate is given by, 

/I i 
d 2 z — — ■ p(z, z* 7 m, p) = G{z = m,z* = to, to, p) (8) 
V Z + TO 

we have solved the Silver Blaze problem if we can show that the spectral density 
satisfies the constraint (JJJ) with G(z = m,z* = to, to, p) independent of p. At p = 0, 
the discontinuity of the chiral condensate is due to an accumulation of eigenvalues on 
the imaginary axis (Banks-Casher relation) . At p ^ the QCD Dirac operator does 
not have any hermiticity properties. The eigenvalues are scattered in the complex 
plane and, because of the phase of the fermion determinant, the eigenvalue density 
is complex. The solution of the Silver Blaze problem will provide us with different 
relation between the discontinuity of the chiral condensate and the spectral density. 
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4. The Dirac Spectrum at Nonzero Chemical Potential 

The spectrum of the quenched Dirac operator is well understood ^ 1 ^ l 1 7 | lS | 22 | l JMil 

In the mean field limit limit at fixed fj, -C Aqcd the eigenvalues are scattered 
uniformly in the strip |Re(z)| < 2/x 2 F 2 /E. The mean field limit of the quenched 
resolvent is therefore given by 

GQ(^*,M) = ^? for |Re(z)|<2 / ?F 2 /S, (9) 

G Q (z,z*,fi) = S for |Re( 2 )| > 2^ 2 i^ 2 /E. (10) 

In the unquenched case we do not have a simple expression for the resolvent except 
that it has to satisfy the constraint J7J at zero temperature. 

The spectral density of the Dirac operator can be obtained fronBSED 

p(z,z*,m,n) = lim — d z d z * log Z n (z, z* , to, fi) (II) 

with generating function Z n (z, z*, (i) defined bjEIEll 

Z n (z, z*,m,fj.) = (det(D + (ijo + m)det™(D + ^70 + z)det n (.D' + (ijo + z*)) 

= (dct(£> + /170 + m)det"(L> + ^ 7o + z)det"(L> - ^ 7o + z*)). (12) 

In general, the replica limit n — > is not well defined^. However, for m,\z\ 
F 2 /(Y 1 t/V) and \i Aqcd this family of partition functions satisfies the Toda 
lattice hierar chy which allows one to o btain the n — * result fro m a recursion 
relatiorPES 

both in the quenchecP^ and in the unquenched cas d^ l The un- 
emenched spectral density, which was first obtaineo from an equivalent random 
matrix model, is a strongly oscillating complex function (see Fig. 3). The amplitude 
of the oscillations grows exponentially with fj, 2 F 2 V while the period goes as 1/Y,V, 
It can be showrJ=2l that the resolvent, obtained by integrating the oscillating spec- 
tral density, satisfies the constraint that G(z — m, z* — to, to, /x) is independent of 
/J.. The discontinuity of the chiral condensate is /i-independent with contributions 
from the entire oscillating regional This solves the Silver Blaze problem. 

The spectral density can also be evaluated in the mean field limit. Then the 
resolvent does not depend on the number of 

replieaJSUMi 

and the replica limit 

can be taken trivially. Assuming that this procedure gives the correct result we 
simply consider n — 1. The mean field generating function (|12f) for n = 1 and real 
z was analyzed before^! by means of a chiral Lagrangian. Interpreting to as the 
strange quark mass and z = z* as the light quark masses, the valence "pion" and 
"kaon" masses are given by, 

2 2xT, 2 (s + to)E 
m V7r = —j- and m vK = — with x = Re(z) (13) 

respectively. With this identification, the isospin chemical potential is equal to fii = 
1\x and the strangeness chemical potential \i s = —fj,. The phase diagram in the x — fi 2 
plane for Im(z) = is shown in the left panel of Fig. 3. We observe three different 



Chiral Symmetry Breaking 5 




Re[ P ]/Y?V 2 xT,V 

Fig. 3. Phase diagram of Dirac spectrum (left) and real part of the spectral density for full QCD 
at nonzero chemical potential (rig . In the left figure full lines denote a second order transition 
and dashed curves a first order transition. The dotted curve (x = §At 2 F 2 /£ - ^) represents the 
right edge of the oscillation region given in the right figure. 

phases, a normal phase (N), a pion condensation phase (jr) and a kaon condensation 
phase (K) bordered by the dashed curve and the solid curve x = 4/i 2 F 2 /E — m. At 
fixed chemical potential, a transition from the 7r-phase to the X-phase takes place 
when m v n = m v k ■ Evaluating the mean field resolvent in each of the phases we 
find 

= E, Gk — E, Gtt — ^ 2 ^ 2 . (14) 

The resulting spectral density is zero in the ./V-phase and a constant in the 7r-phase 
in agreement with the thermodynamic limit of the exact spectral density. However, 
the mean field spectral density in the K-ph&se is zero rather than oscillatory. The 
discontinuity of the mean field resolvent across the interface of the 7r-phase and the 
if-phase gives a delta function in the spectral density. As the quark mass approaches 
zero, the interface approaches zero resulting in a discontinuity of the chiral conden- 
sate due to an accumulation of eigenvalues. This solves the Silver Blaze problem for 
the mean field spectral density for Nf = 1 obtained from n = 1. 

5. Discussion 

We have described two different mechanisms to obtain a discontinuity in the chiral 
condensate at nonzero chemical potential. First, in case of the low energy limit of 
the full QCD partition function, the discontinuity arises from the oscillating part 
of the spectral density. Second, using a replica mean field result for the spectral 
density we find that the discontinuity moves with m and reaches the imaginary axis 
for m = 0. At this moment it is not yet clear whether or how these two pictures 
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are related. One important observation is t hat the oscillating region of the spectral 
density is absent for fi < We expecJ^that the sign problem is less severe 

this in this reg ion of the phase diagram which could open the door to realistic lattice 
simulations^! 
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